Coarse-Graining and Renormalization of Atomistic 
Binding Relations and Universal Macroscopic 

Cohesive Behavior 



O. Nguyen and M. Ortiz 
Graduate Aeronautical Laboratories 
California Institute of Technology 
Pasadena CA 91125, USA 



Abstract 

We present two approaches for coarse-graining interplanar potentials and 
determining the corresponding macroscopic cohesive laws based on energy re- 
laxation and the renormalization group. We analyze the cohesive behavior of 
a large — but finite — number of interatomic planes and find that the macro- 
scopic cohesive law adopts a universal asymptotic form. The universal form of 
the macroscopic cohesive law is an attractive fixed point of a suitably-defined 
renormalization-group transformation. 

1 Introduction 

Cohesive theories of fracture are predicated on a direct description of the phys- 
ical processes which lead to separation and the eventual formation of a free 
surface. The development of cohesive theories rests on a detailed physical un- 
derstanding of the operative fracture mechanisms, which are often complex and 
cut across multiple lengthscales, especially where ductile fracture is concerned. 
Cleavage fracture, by way of contrast, entails the simple separation of atomic 
planes and is, therefore, governed by interplanar potentials which are amenable 
to an effective first-principles atomistic characterization. For instance, Jarvis et 
al. H have recently calculated the cohesive behavior of (111) planes in fee alu- 
minum, and of AI2O3 cleavage planes, using GGA density functional theory; 
and Park and Kaxiras @ have carried out ab-initio simulations of hydrogen 
embrittlement in aluminum and calculated generalized stacking-fault energies 
as a function of interplanar separation and sliding. 

First-principles interplanar potentials are characterized by peak stresses of 
the order of the theoretical strength of the crystal. In addition, the crystal loses 
its bearing capacity after an interplanar separation of only a few angstroms. 
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Moreover, the integration of first-principles interplanar potentials into engineer- 
ing calculations necessitates full atomistic resolution in the vicinity of the crack 
tip, which is often unfeasible or impractical. This disconnect between atomistic 
and engineering descriptions begs a number of fundamental question, to wit: 
What is the proper way to coarse-grain a cohesive description?, and: What is 
the macroscopic form of the cohesive law after coarse-graining? 

In this paper we address these issues by investigating the cooperative be- 
havior of a large number of interatomic planes forming a cohesive layer. We 
employ two main approaches in this investigation: relaxation and the renor- 
malization group. Relaxation or weak convergence methods are concerned with 
the determination of the macroscopic behavior of materials characterized by 
a non-convex energy function. These materials often develop fine microstruc- 
ture in response to imposed deformations. Truskinovsky et al. |3|, |8|, Q, and 
Braides et al. ||, have pioneered the application of these methods to fracture. 
However, the full relaxation of a cohesive potential yields the trivial result that 
the effective cohesive potential is identically equal to zero. The chief difference 
between the analysis pursued here and full relaxation is that, at zero tempera- 
ture, we seek energy minimizers of large — but finite — collections of interatomic 
planes. In this limit we find that, for a broad class of interplanar potentials, 
the macroscopic cohesive law adopts a universal form asymptotically. 

We show that this universality of the macroscopic cohesive behavior is 
amenable to a renormalization-group interpretation. The normalization group 
which coarse-grains the cohesive behavior is somewhat nonstandard and has to 
be crafted carefully, e. g., so as to preserve the surface energy and the elasticity 
of the lattice. The universal form of the macroscopic cohesive law is precisely 
an attractive fixed point of the renormalization-group transformation. 



2 Problem formulation 

We consider a macroscopic cohesive crack opening symmetrically (mode I) and 
undergoing quasistatic growth. We denote by d the interplanar distance, 5 
the opening displacement across an interatomic plane, and t the corresponding 
cohesive traction. These latter variables are presumed related by a known 
cohesive or binding law t(5), which derives from an interplanar potential (f>(5) 
through the relation 

W = (!) 

Here and subsequently, a prime denotes differentiation of a function of a single 
variable. For simplicity, we shall assume throughout that the atomistic binding 
law t(5) rises monotonically from zero at 5 = to a peak value a c at 5 = 5 C , 
and subsequently decreases monotonically to zero, Fig. |l]. Correspondingly, the 
cohesive potential (f)(5) is convex in the interval < 5 < 5 C , has an inflection 
point at 5 = 5 C , is concave for 5 > 5 C and asymptotes to twice the surface 
energy, 27, as 5 — ► 00. In addition, we shall assume that (f)(5) is smooth and 
analytic at 5 = 0, with Taylor expansion: 

cp~^5 2 + o(5 2 ) (2) 
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for some constant C. The value of C can be readily deduced from the elastic 




Figure 1: Interplanar potential and corresponding cohesive or binding law. 



moduli Cijki of the crystal. To this end, let m be the unit normal to the 
plane of the crack, and apply a small and uniform opening displacement to all 
interatomic planes. Evidently, the energy per unit volume of the crystal follows 
from the cohesive potential as (C/2d)5 2 asymptotically as 5 — > 0. On the other 
hand, the strain tensor of the crystal is = (5/d)mimj, and the corresponding 
energy is (5 2 /2d 2 )cijkimimjmkmi. Equating both energies yields the identity: 



n 1 



(3) 



Next imagine that the atomistic description is coarse-grained, e. g., by the 



quasi-continuum method [10|, or by a passage to the continuum limit, or by 
some other suitable means. Let d denote the spatial resolution of the coarse- 
grained description. For instance, in quasicontinuum or in engineering finite- 
element simulations d measures the local element size. The corresponding ef- 
fective cohesive law may be obtained by analyzing the behavior of a cohesive 
layer of thickness d and containing N = d/d atomic planes. The cohesive layer 
is taken through a total opening displacement 5 resulting in a macroscopic 
traction t. The chief objective of the analyses that follow is to determine the 
macroscopic cohesive law t(5) in the limit of N large but finite. Equivalently, 
we may seek to determine the asymptotic form of the macroscopic cohesive 
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potential (f)(8) such that 

t(6) = 4>'(6) (4) 

in the same limit. 

3 Universal asymptotic form of the macro- 
scopic cohesive law at zero temperature 

At zero temperature, the crystal deforms so as to minimize its total energy. The 
governing principle is, therefore, energy minimization. Let 5j > 0, i = 1, . . . , N 
be the opening displacements of the interatomic planes in the cohesive layer. 
Then, the total energy of the cohesive layer is: 

N 

£ tot = J>(^) (5) 

8=1 

Let now 5 be the macroscopic opening displacement. Then, the effective or 
macroscopic energy of the cohesive layer follows from the constrained mini- 
mization problem: 

N 

W = r^H (6) 
N 

5 = J> (7) 

i=l 

In conjunction with the kinematic constraints (0) the stationarity of E tot de- 
mands: 

t(5i) = i(S) = constant, i = 1, . . . , N (8) 

Thus, at equilibrium, all interplanar tractions must be equal to the macroscopic 
traction. 

We shall classify the possible states of an interatomic plane into two cat- 
egories or variants, according as to whether the opening displacement 5 is in 
the range < 5 < 5 C , or in the range 5 > 5 C . We shall designate variants of 
the first kind as coherent, and variants of the second kind as decohered. We 
may further classify the states of a cohesive layer by the number N\ of coher- 
ent planes, or, equivalently, the number ./V2 of decohered planes, it contains. 
Since the function t(5) is one-to-one over the interval [0, 5 C ), eq. (||) demands 
that the opening displacements of all coherent planes be equal at equilibrium. 
Likewise, the opening displacements of all decohered planes must be identical 
at equilibrium. Under these conditions the macroscopic cohesive energy follows 
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from the minimization problem: 

m = m mwj+Nrifr)} (9) 

{{oi,02),{Nl,N2)i 

6 = N 1 5 l + N 2 8 2 (10) 

N = N± + N 2 (11) 

< Si < 5 C (12) 

S 2 > 6c (13) 

where 5i and 6 2 are the opening displacements in the coherent and decohered 
planes, respectively. In addition, the equilibrium equations (||) reduce to: 

= t(5 2 ) = i(S) (14) 

These relations are depicted geometrically in Fig. [|. 

Next we proceed to determine the minimum energy states of a cohesive layer 
by analyzing the cases N 2 = 0, 1,2, . . . in turn. We begin by considering the 
case in which all planes are coherent, corresponding to N\ = N and N 2 = 0. 
Then, the kinematic constraint (|^) gives 5\ = 5/N = 5. Evidently, in the limit 
of N — > 00 5 tends to zero and, in view of eqs. @ and @, we obtain: 

^)|j Va =o~? 52 ' asAr ^°° ( 15 ) 

where 

is an effective cohesive-layer stiffness. 

Consider next the case of one decohered plane, N± = N — 1 and N 2 = 1, 
whence (|l0|) becomes: 

(N - 1)5! + 6 2 = 6 (17) 

Solving for <5i gives 5± = (5 — 5 2 )/(N — 1). In addition, since 5i and 8 are 
required to be nonnegative, it follows that 5\ < 5/(N — 1), and thus 5\ — > 
as N —* 00. From this limit it additionally follows that 5 2 — > 5 in the same 
limit. Suppose now that <5 3> <5 C and, hence, S 2 3> 5 C . Under these conditions, 
4>{5 2 ) ~ 27 and (||) reduces to: 

^)|at 2= i~ 2 7> asiV^oo (18) 

Since interactions beyond nearest neighbors are not taken into account, an 
altogether identical analysis gives 

$(8)\ N] 2=k ~ 2^7, asiV^oo (19) 

for the case of A; decohered planes. 

The macroscopic cohesive energy may now be expressed as 

j>(5) = min 0(6) I Ar , (20) 
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However, it follows from ([11]) that, asymptotically as N — > oo, multiple deco- 
hered planes are not energetically possible at zero temperature, and only the 
cases k = and k = 1 need be considered in (|20|) . Therefore, the effective 
cohesive potential is the lower envelop of the energies (15) and (18), namely, 



■ r C Sm f (<?/2)J 2 , if5<<5 c , ni , 
m=mm{-fi 1 2n,} = [ /} otherwige (21) 



where 



*= 2 Vl= 2 vf (22) 

is a macroscopic critical opening displacement for the nucleation of a single 
decohered plane. The corresponding macroscopic cohesive law is 

It is interesting to note that the peak macroscopic traction is: 

a c = C6 C = 2^ = 2^- (24) 

We conclude the analysis by verifying that, in the decohered regime, 5 > 5 C ~ 
y/~N, and hence 5 3> 5 C for sufficiently large N, as supposed. 

These functions are shown in Fig. |2[ The macroscopic cohesive potential is 
initially quadratic and subsequently constant following the attainment of the 
critical macroscopic opening displacement. Remarkably, the macroscopic criti- 
cal opening displacement and peak traction scale as: 5 C ~ \/iV and c c ~ 1/ \/iV, 
respectively. Thus, for large JV, it follows that the macroscopic cohesive law 
entails much lower tractions, occurring at much larger opening displacements, 
than the atomistic binding law. In effect, the passage from the atomistic to the 
macroscopic scales is accompanied by an expansion of the opening displacement 
axis and a simultaneous compression of the traction axis. By constrast, the 
macroscopic fracture energy, or critical energy-release rate, 4>(oo) remains in- 
variant under the transformation and is equal to the atomistic value <f>(oo) = 27. 
It is also remarkable that, for the class of binding laws under consideration, the 
asymptotic form (|2l|) of the macroscopic cohesive law is universal, i. e., inde- 
pendent of the atomistic binding law. Evidently, the parameters which define 
the macroscopic cohesive law quantitatively, e. g., the surface energy 7 and the 
modulus C, are material specific. 

As a simple illustrative example we consider the universal binding energy 
relation (UBER) || defined by the interplanar potential: 

<P{5) = 2 7 - CS C (5 + 5 c )e- 5 ^ (25) 

This function falls within the class of potentials considered in the foregoing. 
We choose as material constants: C = 3.54 J/m 2 /A 2 , S c = 0.66 A, which are 
representative of aluminum. The macroscopic cohesive laws resulting from a 
direct numerical minimization of the energy (w for different values of iV are 
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Figure 2: Universal asymptotic form of the macroscopic cohesive law at zero temper- 
ature. 



shown in Fig. The universal asymptotic form of the macroscopic cohesive 
law is compared in Fig. ||| against the corresponding numerical results. The 
convergence of the macroscopic cohesive law towards the universal asymptotic 
form as the number of planes in the cohesive layer increases is clearly evident 
in this figure. 

4 Effect of finite temperature 

At finite temperature, entropic effects make it feasible for the cohesive layer 
to decohere on multiple planes. In order to assess this effect simply we recall 
that, asympotically, the energy of a cohesive layer with no decohered planes is 
(C/2) 5 2 , and that the energy of a cohesive layer containing k decohered planes 
is 2Ar/. Within this approximation, the partition function of an area a 2 of layer 
is, therefore, 

oo 

Z(S, T) = e-^^~ s2a2 + J2 ^ 2kW (26) 

k=l 

where (3 = 1/kT, and k is Boltzmann's constant. In order to count states 
properly, we identify a with the lattice parameter of the crystal. Physically, 
this is tantamount to allowing for decohered areas of a size commensurate with 
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Figure 3: Numerically computed macroscopic traction vs opening displacement rela- 
tion for an increasing number of atomic planes in the cohesive layer. 



the lattice parameter. The sum in (|26| ) defines a geometric series which may 
be evaluated readily, with the result: 

- -„ „ -/327a 2 

zm = e ^c /2)S ^ + _^__ (27) 

The free energy density per unit area of the layer now follows as 

F(6,T) = ~~]ogZ(5,T) (28) 
o z p 

whereas the resulting effective cohesive law is: 

dF - 

i(5,T) = ^(S,T) (29) 

The effect of temperature on the macroscopic cohesive potential for alu- 
minum, endowed with an interplanar potential of the UBER type, is shown 
in Fig. The lattice parameter is taken to be a = 4.05 A. As expected, 
the Helmholtz free energy rises above the zero-temperature in the amount TS, 
where S is the configurational entropy of the layer, Fig. ||a. The correspond- 
ing effect on the macroscopic cohesive law is to smooth out the decohesion 
transition, Fig. |B|b. 

5 Renormalization Group interpretation 

The renormalization group (RG) (see, e. g., ||) provides a natural framework 
for the understanding of universality, i. e., the phenomenon that large classes of 
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(c) N = 1000 
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Figure 4: Comparison between numerically computed macroscopic cohesive law and 
universal asymptotic form for an increasing number of atomic planes in the cohesive 
layer. 
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systems with unrelated Hamiltonians may nevertheless exhibit identical ther- 
modynamic behavior near critical points. It is therefore not entirely unexpected 
that the main result of Section ||, namely, that the limiting form of the potential 
of a cohesive layer is universal for a broad class of interplanar potentials, can 
be given a compelling interpretation within the RG framework. 

For simplicity, we confine our attention to the zero temperature case. We 
proceed to construct an RG transformation R such that the sequence 

4> n+ i = R(f> n , n = 0, 1, . . . (30) 

with (fio (6) = 4>(5), yields, by recourse to an appropriate scaling, the large N 
asymptotic form of the macroscopic cohesive law in the limit. As in the pre- 
ceding section, the interplanar potentials contemplated here, and to which the 
transformation R is applied, are continuous, monotonically increasing functions 
<f> : [0, oo) — > [0, 27] and analytic at the origin. For simplicity, we additionally 
restrict attention to functions (f)(6) possessing a single inflection point, so that 
4>'(6) has a single maximum. 

We construct R by the usual combination of decimation and scaling. The 
decimation step concerns a layer consisting of two interatomic planes with open- 
ing displacements 5\ and £2 and total opening displacement 8. The correspond- 
ing effective energy follows from the minimization problem: 

m = inf {^1) + ^)} (31) 
{01,02} 

6 = 5i + 6 2 (32) 

Next, we proceed to rescale (j)(5) in such a way that the sequence defined by 
d30| ) has a well-defined fixed point. Since the transformation must preserve the 
relation (f)(oo) = 27, it is clear that we are allowed to rescale the independent 
variable 6 only. Thus, we set: 

(R(f>)(6) = 4>(\5) (33) 

For very small 8 the interplanar potential (f> is essentially quadratic and reflects 
the elasticity of the lattice. Thus, in order for the transformation to preserve 
the elasticity of the lattice it must leave parabolic functions invariant. For <p = 
(C/2)5 2 it follows that <5i = 5 2 = 8/2 and 4> = 2{C/2)(6/2) 2 = (C /2)(5 / s/2) 2 . 
Finally, (R(f))(8) = (C/2)(A<5/ v / 2) 2 , whence it follows that A = a/2 for to 
remain invariant under the transformation. The complete RG transformation 
is, therefore, 



4>(6) = inf U(6 1 ) + <P(6 2 )} (34) 
{01,02} 

6 = 5 1 + 5 2 (35) 

(R4)(S) = 4>(V26) (36) 



Taking 6\ = and 6 2 = 6 in (^J) immediately shows that the unsealed energy 
(f)(6) is bounded above by the original function (f)(6). 

It is clear from definition ( |34[|3"6| ) that the RG transformation leaves the 
specific fracture energy (f)(oo) invariant and equal to its initial value 27. Another 
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invariant of the RG transformation is the initial modulus C = <p"(0). Indeed, 
consider the limit of R<j) as 5 — ► 0. Since necessarily 5\ < S and 62 < 5 it follows 
that both 5\ — > and 82 — > in this limit. (R(f>)(8) may therefore be computed 
by replacing </>(<5) by its Taylor expansion about the origin, namely, (C/2)8 2 , 
with C = 0"(O). But parabolic functions are invariant under R and, hence, so 
is C. 

The RG transformation R (|3^ - ^6|) preserves the monotonicity of </>(<5). In 
order to see this, consider a pair of opening displacements 8 and 5' = X8, with 
A < 1. Let (f>(5) = (f)(5i) + ^(62) for some pair of opening displacements 81 
and 82 satisfying constraint (|35l). The opening displacements 8[ = X5\ and 
8' 2 = X82 then satisfy the similar constraint: 8' = 8[ + 8' 2 - It therefore follows 
that 4>(5') < 4>{X5i) + <p(\S 2 ) < <t>{8\) + <fi(5 2 ) = 4>(8). An application of the 
rescaling ( |36| ) to both sides of this inequality finally proves the assertion. 

It is easy to show that the function: 

(«5) =min{^ 2 ,2 7 } (37) 



c 



is a fixed point of i?. To this end, we may distinguish the cases: a) 5± < 5, 
and 62 < 5 C ; b) 5\ < 5 C and 82 > <5 C , or, equivalently, 62 < <5 C and <5i > <5 C , and 
c) Si > S c and 52 > <5c- Case (a) requires that <5 < 25 c and gives an unsealed 
energy: (f)(5) = (C /2)(d/\/2) 2 . Case (b) requires that 5 > 5 C . The optimal 
unsealed energy is obtained by taking Si = and 62 = 5, with the result: 
(f)(8) = 27. Case (c) results in the unsealed energy: (f)(5) = 47. The function 
j37| ) is recovered by taking the minimum of the unsealed energies resulting from 
cases (a), (b) and (c) and applying the scaling (p6| ) to the result. 

A key question is whether the fixed point ( |37D is attractive. We have in- 
vestigated this question numerically for the particular example of the UBER 
binding law, (|25|) . Fig. || shows the evolution of (f> n with increasing n. It is clear 
from the figure that, at least for the example under consideration, the flow of 



functions (f> n (S) does indeed converge strongly to the fixed point (37). The re- 
lation to the asymptotic cohesive law (^Tj) is as follows. We may regard (f> n (5) 
as the result of decimating and rescaling n times a cohesive layer containing 
N = 2 n planes. The total opening displacement of the layer is obtained by 
undoing all the rescalings, with the result: S = (y/2) n 5 = y/NS. For large N, 
<t>n(S) ~ <j>oo($) and one has 



</>{S) ~ <f>oo(S/VN), as N -> 00 (38) 



which is identical to (|2l[). 

Eq. (p^) establishes a connection between the renormalization group, specifi- 
cally as generated by transformation ( ^4[j36| ), and the large- ./V asymptotic form 
of the cohesive potential determined directly in the preceding section. It is 
interesting to note that the RG transformations which pertain to the renormal- 
ization of interplanar potentials are markedly different from those which arise 
in the calculation of bulk thermodynamic properties. In this latter context, the 
appropriate scaling is related to the volume of the sample and is designed so as 
to result in well-defined extensive fields and intensive variables. In the present 
context, the energy densities under consideration are defined per unit area, and 
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the limit of interest is the total energy of the cohesive layer per unit surface 
area, as opposed to the energy per unit volume. In addition, the indepen- 
dent variable of interest is the total opening displacement across the cohesive 
layer, as opposed to its transverse strain. These peculiarities account for the 
non-standard character of the renormalization group defined in the foregoing. 

6 Summary and conclusions 

We have presented two approaches for coarse-graining interplanar potentials 
and determining the corresponding macroscopic cohesive laws based on energy 
relaxation and the renormalization group. We have analyzed the cohesive be- 
havior of a large — but finite — number N of interatomic planes and found that 
the macroscopic cohesive law adopts a universal asymptotic form in the limit of 
large N. We have also found that this asymptotic form of the macroscopic cohe- 
sive law is an attractive fixed point of a suitably-defined renormalization-group 
transformation. 

The universal asymptotic form of cohesive law is particularly simple: the 
traction rises linearly from zero to a peak stress a c at a critical opening dis- 
placement 6 C , and subsequently drops to zero. The scaling of the peak stress 
and critical opening displacement is 5 C ~ 1/y/N and 5 C ~ yN. Thus, coarse- 
graining is accompanied by an attendant reduction (increase) in the cohesive 
traction (opening displacement) range, while at the same time preserving the 
surface of specific fracture energy of the crystal. 

It is interesting to note that the size of a cohesive zone at the tip of a 
crack in an elastic crystal scales as I ~ l/°c 0- Upon coarse-graining, the 
cohesive zone size increases to / ~ l/a c 2 , which gives I ~ Nl. This scaling 
preserves the ratio 1/1 = d/d, which shows that coarse-graining has the effect of 
expanding the cohesive-zone size to within the resolution of the coarse-grained 
description. In particular, it eliminates the onerous need to resolve the atomic 
scale in simulations. 

It is also interesting to note that the universal form of the macroscopic 
cohesive potential is completely determined by the constants C = 4>"(0) and 
4>(oo) = 27. This greatly reduces the scope of the first-principles calculations 
required to identify the macroscopic cohesive behavior of specific materials, 
which can be limited to the calculation of elastic moduli, lattice constants and 
surface energies. 

Finally, we close by suggesting possible extensions of the theory. The anal- 
ysis presented in the foregoing has been restricted to symmetric (mode I) open- 
ing normal to the atomic planes. A worthwhile extension would be to consider 
interplanar potentials defined in terms of three opening displacements and, 
therefore, capable of describing tension-shear coupling. Another worthwhile 
extension would be to consider interplanar potentials with multiple inflection 
points, which would greatly enlarge the class of materials tractable within the 
theory. 



O. Nguyen and M. Ortiz 



Coarse-Graining and Renormalization. . . 13 



Acknowledgements 

This work has been supported by Brown University's MURI Center for the 
"Design and Testing of Materials by Computation: A Multi-Scale Approach." 
We are grateful to Emily A. Carter, Emily A. A. Jarvis and R. L. Hayes for 
many useful discussions and suggestions, and for making available to us their 
research results prior to publication. 



References 

[1] G. I. Barrenblatt. The mathematical theory of equilibrium of cracks in 
brittle fracture. Advances in Applied Mechanics, 7:55-129, 1962. 

[2] A. Braides, G. Dal Maso, and A. Garroni. Variational formulation of soft- 
ening phenomena in fracture mechanics, the one-dimensional case. Archive 
for Rational Mechanics and Analysis, 146(l):23-58, 1999. 

[3] G. del Piero and L. Truskinovsky. A one-dimensional model for localized 
and distributed failure. Journal de Physique IV, 8(P8):95-102, 1998. 

[4] G. Del Piero and L. Truskinovsky. Macro- and micro-cracking in one- 
dimensional elasticity. International Journal of Solids and Structures, 
38(6-7): 1135-1148, 2001. 

[5] N. Goldcnfcld. Lectures on Phase Transitions and the Renormalization 
Group. Addison- Wesley, 1992. 

[6] E.A.A. Jarvis, R.L. Hayes, and E.A. Carter. Effects of oxida- 
tion on the nanoscale mechanisms of crack formation in aluminum. 
Chem. Phys. Chem., 2(l):55-59, 2001. 

[7] I. Part and E. Kaxiras. Ab-initio studies of hydrogen embrittlement in 
aluminum, to appear, 2001. 

[8] G. Puglisi and L. Truskinovsky. Mechanics of a discrete chain with bi- 
stable elements. Journal of the Mechanics and Physics of Solids, 48(1):1- 
27, 2000. 

[9] J.H. Rose, J.R. Smith, and J. Ferrante. Universal features of bonding in 
metals. Physical Review B-Condensed Matter, 28(4):1835-1845, 1983. 

[10] E.B. Tadmor, M. Ortiz, and R. Phillips. Quasicontinuum analysis of de- 
fects in solids. Phil. Mag. A, 73(6):1529-1563, 1996. 



O. Nguyen and M. Ortiz 



Coarse-Graining and Renormalization. 



14 



2 5 




Macroscopic separation [Angstroms] 

(a) Helmholtz free energy. 




(b) Macroscopic cohesive traction. 



Figure 5: Influence of the temperature on the effective behavior of the chain (TV = 
100). 
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